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Abstract 



We study the transport of an overdamped particle adiabatically driven by 
an asymmetric potential which is periodic in both space and time. We develop 
an adiabatic perturbation theory after transforming the Fokker-Planck equa- 
tion into a time-dependent hermitian problem, and reveal an analogy with 
quantum adiabatic particle transport. An analytical expression is obtained 
for the ensemble average of the particle velocity in terms of the Berry phase 
of the Bloch states. Its time average is shown to be quantized as a Chern 
number in the deterministic or tight-binding limit, with exponentially small 
corrections. In the opposite limit, where the thermal energy dominates the 
ratchet potential, a formula for the average velocity is also obtained, showing 
a second order dependence on the potential. 

PACS numbers: 05.40.-a, 02.40.-k, 72.10.-d, 87.10.+e 
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The so-called ratchet models have attracted much interest in recent years (T|-§[. In a 
typical ratchet model, a driving force with time-correlated fluctuation, in addition to a 
random force, is exerted on an overdamped particle. A macroscopic flow develops if the 
parity symmetry is broken in the driving potential. This behavior is in striking contrast to 
the equilibrium situation in which the second law of thermodynamics forbids a non- vanishing 
macroscopic flow, irrespective of the microscopic details ||. Here we consider an adiabatic 
ratchet potential with both spatial and temporal periodicities, which has received some 
attention recently 0J^]. We take an analytical approach by transforming the Fokker-Planck 
equation into a time- dependent hermitian problem. The ensemble average of the particle 
velocity is calculated by using the standard adiabatic perturbation theory, whose validity 
condition is also explicitly specified for the present system. The particle transport in a 
time cycle is found in terms of the Berry phase of the Bloch waves. In the deterministic 
or tight-binding limit, in which the potential dominates the thermal energy, the particle 
transport is shown to be quantized as a Chern number [^], consistent with the previous 
numerical result ||, and in analogy with quantum adiabatic particle transport first studied 
by Thouless and Niu P-jlOfl. The correction to the quantization is found to be exponentially 
small based on an analysis using the Wannier functions. We also elaborate the opposite 
limiting case, where the thermal energy dominates the potential. 

Consider an overdamped particle moving in one-dimension, subject to a driving force 
F(x,i) = —d x V(x,t) and a random force u(t). It obeys a Langevin equation, 

^ = -F(x,t) + u(t), 
at jj, 

where 7 is the friction constant and p is the mass. We assume that v(t) is of zero average and 
5-correlated, i.e. (y(t)v(t')) = 2(^y ksT / p)5(t — if) , where ks is the Boltzmann constant, T is 
the temperature, and (■ ■ •) denotes the ensemble average. The potential V(x,t) is periodic 
in both space and time, with periods a and r respectively. Furthermore, the potential is 
assumed to be asymmetric, i.e. V(x) 7^ V( shown in Fig. 1. 

The corresponding Fokker-Planck equation for the probability density p(x, t) is 

-d tP (x,t) = DOp{x,t), (1) 

where D = ksT '/ '/ry is the diffusion constant, and O = — d^+d x - F/kaT. The Fokker-Planck 
equation ([I]) can also be written as a continuity equation d t p+d x j = 0, where the probability 
current is j(x,t) = Djp(x,t), with j = —d x + F/k B T. 



The operator O is non-hermitian, but there is a similarity transformation [TJ]] leading to 
a hermitian operator, to be called the effective hamiltonian, 

H = e r Oe~ r =p 2 + U, (2) 

where T = V/2ksT is the potential normalized by the thermal energy, U = {d x T) 2 — d 2 x Y 
is the effective potential, and p = —id x . Note that 7i has the dimension of inverse length 
squared. Since 7i is hermitian, its eigenfunctions comprise a complete orthonormal basis. 
These are the Bloch waves, ip n k, owing to the periodicity of the potential, where n is the 
discrete band index, and k is the Bloch vector. It is also interesting and useful to note that 
the ground state eigenvalue is zero, and the wave function can be explicitly found as 

/ 1 -r 
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where Z = J e 2T dx is a normalization factor ensuring (V'olV'o) = / i/Jo(x)*ipo(x)dx = 1. 
The Fokker-Planck equation is then transformed into 

- d t ip(x, t) = (DH - d t T - d t In \[Z)^{x, t), (3) 

where ip = pe r \fZ is the transformed variable for the probability density, and the factor 
\/Z is added for convenience. If ip = ipQ, then the probability density is the stationary 
Boltzmann distribution p = e~ 2T /Z, giving j = 0. If V were time-independent, the above 
equation would reduce to — d t ip = DTCip, whose ground state solution ip = ip corresponds 
to the equilibrium state. 

The wave function ip is generally a superposition of different Bloch states. However, 
only those with k=0 are important if we are only interested in the ensemble average of the 
particle velocity (%). This is equal to the spatial integral of the probability current, 

J{t) = D J jpdx = D J ip (-d x - d x Y)ij)dx = -2D J i/) d x i/)(x,t)dx, (4) 

where, in the expansion of ip(x,t), only if} n k{sc,t) with k = contribute. This result does 
not require the necessity of adiabaticity. 

We now apply the adiabatic perturbation theory [0 , supposing that the ratchet poten- 
tial changes slowly. We first make the expansion 

1>(x,t) =Y, c n(t)Mx,t)e~ Df ° m ' )dt ', (5) 

n 

where k = is assumed in light of the discussions above. For clarity, we have retained 
i?o although it is zero. The coefficient for the ground state is constrained to Co(i) = 1 as 
imposed by the normalization of the density p(x,t). The coefficients for the higher bands 
satisfy 

D{E - E n )c n + ]T cv(V> n |<%r + d t In ^[Z - <%|Vv) = d t c n . (6) 

n> 

The terms in the summation are small if the potential changes slowly and if we choose the 
Bloch waves with k = to be real such that < ip n \dt\?pn >— 0. We also assume that the 
system is initially in the ground state (stationary state). Thus we can ignore all but the 
n' = term in the summation, obtaining, to first order in the rate of change of the potential, 

° nMt) D[E (t) - E n (t)Y [() 

where we have used the fact that (dtT + d t In y/Z)ip = —d t ipo- The smallness of c n for n > 
means that the system remains close to quasi-static equilibrium with the instantaneous 
ratchet potential. If the system has initial excitations in the higher bands, one can show 
that they decay away exponentially at a rate of D(E n — Eq) and become negligible after a 
transient period. 

Quantitatively, the adiabatic condition is satisfied if the potential changes at a rate 
much smaller than D times the eigenvalue gap, - D(Ei — E ). We now make an order of 
magnitude estimate of this gap. We assume that the amplitude of variation in the force is 
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Fq, yielding the variation for T as To = -Foa/2/csT. In the kinetic regime where To << 1, the 
gap at k = is about (27r/a) 2 . The adiabatic condition is just r ^> t d , where r D = a 2 / D is 
the diffusion time over one wavelength of the ratchet. In the tight-binding or deterministic 
regime, where T >> 1, the Bloch bands may be viewed as derived from the levels in the deep 
potential wells. A peculiar thing to notice is that the effective potential U typically contains 
a double well structure in each wavelength of the ratchet, even if the ratchet potential has 
one well in each period. The reason is that the dominant term (d x T) 2 in U has only half 
of wavelength of the ratchet, and thus must contain two wells in each wavelength of the 
ratchet. These two wells are made inequivalent by the weaker term, —d 2 T, which has the 
full periodicity of the ratchet. The band gap can thus be estimated by the level difference 
between the two wells, and is given by the amplitude of variation of the weaker term To/a 2 . 
The adiabatic condition in this case is r ^ td/Tq, which is more easily satisfied than the 
kinetic regime because of the large T in the tight-binding regime. 

The average velocity induced by the adiabatic movement of the ratchet can then be 
obtained from (f|), @ and (0) as 

This current is independent of the diffusion constant D, because all the eigenstates and 
eigenvalues of the effective Hamiltonian are independent of D. This signifies that the adia- 
batic flow is a geometric effect: it only depends on how the ratchet potential evolves in time, 
and is independent of the mobility of the particles. The only role of the diffusion constant 
is to set the time scale for the validity of the adiabatic approximation. The reader is also 
reminded that the distinction between the tight-binding and kinetic regimes, in which the 
adiabatic flow will be shown to have very different behaviors, also has nothing to do with 
D. 

In the following, we make a further connection to the Berry phase, and show when the 
adiabatic current may be quantized. We can write J(t) = J{k = 0, t), with 



J(k,t) = -2J2 



Eot — E nk Eok — E nk 



(9) 



This can be rewritten in terms of the periodic amplitude u n k{x,t) of the Bloch waves, and 
turned into a form involving the lowest band only 

J(k,t) = ~ l(d k u ok \d t u ok ) - (dtu 0k \d k u 0k )} , (10) 

following the method used in deriving the quantum adiabatic particle transport ||,[|. This 
is known as the Berry curvature of the Bloch state in the parameter space of k and t. The 
time average of J(t) over a time cycle of the ratchet is therefore 

1 rr 1 

J(t) = — dt {(d k u 0k \dtu 0k ) - (d t u ok \d k uok)} k= Q = — d k 9(k) k=0 , (11) 

IT JO " T 

where 9(k) — i j dt(uQ k \d t UQ k ) is the Berry phase for a given k. In obtaining the last ex- 
pression in the above equation, we assumed a phase gauge such that the wave function is 
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periodic in time, which is possible because the eigenstate must come back to itself in a time 
cycle of the potential apart from a phase factor. 

In the tight-binding or deterministic limit, the quantity J{k, t) is insensitive to k, as can 
be seen in (0), which can be be written as 



J(k,t)=2 



(d x ip k\ ^j ldtipok) + (d t tpok\ y-i 1 -, \d x ipok) 
&ok — n boh — n 



(12) 



We may thus replace the k = expression in (11) by an average over k, yielding 



2w/a 



J(t) = —— / dk dt [(d k u ok \d t uok) - (d t uok\dkU k)] ■ (13) 
Zmr Jo Jo 

This is a closed surface integral of the Berry curvature, and can only take quantized values, 
i.e. 



J(t) = N 



a 



where N is an integer called Chern number. 

Corrections to this exact quantization, which arise in our replacing the k=0 expression by 
an average over k, can be shown to be exponentially small in IV Consider first the somewhat 
trivial example of a sliding potential, V{x — vt), which is in general not a ratchet. The 
position and time dependence of the Bloch states must be through the combination x — vt. 
We can thus replace d t by — vd x , turning (8) into an expression related to the "effective 
mass" of a band [|1(J. The result is that the average particle current can be written in the 
form v[l — 2~ 1 d 2 Eok/dk 2 \k=o\. The first term is just the quantization with N — 1, which 
means that the particle follows the potential, while the second term is exponentially small 
in To because the band width is so in the tight-binding limit. 

To establish the exponential smallness of the correction to quantization for the general 
case of the tight-binding regime, we introduce the the Wannier functions <p n {x — la) = 

J dke~ lkla "il) n k, which are known to be exponentially localized and form an orthonormal 
basis for each band. We may rewrite (f[0|) , with k = 0, as 

J(t) = d t (x) + aY,l(Mx,t)\d t Mx ~ M)>> (14) 
i 

where (x) = (<j>o(x, t)\x\4> (x, £)) is the center position of the Wannier function [|K|. In deriv- 
ing the above result, it is useful to recall the inverse relation ip nk (x,i) = ■^jJ2i elkla 4 ) n(x — 
la, t) between the Bloch and Wannier functions, where M is the number of cells in the 
system. Therefore 

Tttj =-((x)\ t=T - (x)\ t=0 ) + -J2 l f ' dt(Mx,t)\dtMx-la,t)). (15) 

T T i JO 

(x) can only change by an integer multiple of the lattice constant in a time period, so the first 
term in ( |15|) is nothing but the quantization Na/r. The second term gives a correction to 
this quantization. Contribution from I = term is automatically zero. In the tight-binding 
limit, these Wannier functions decay exponentially away from the bottom of each well, so the 
leading term in the correction comes from the nearest-neighbor term 1 = 1. Therefore, the 
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correction is proportional to the overlap between the nearest-neighboring Wannier functions, 
which is exponentially small in T . 

On the other hand, beyond the tight-binding or deterministic limit, the quantization is 
destroyed. We now consider the opposite limit, the kinetic regime where T << 1. We 
can make a perturbative expansion the Bloch function in terms of the plane waves, tpok = 
E m Cm(% i(fc+m2l/a)l , where c (k) = 1 and c m = U m /[k 2 - (k + 2nm/a] 2 ) for m ^ 0, with 
U m = J U(x,t)e~ im27TX ^ a dx being the Fourier coefficient of the effective potential. Because 
we are interested in the neighborhood of k — 0, no degeneracy needs to be considered. The 
average velocity is calculated as J{t) = ^ 27r / a )5 (U-id t Ui — U\d t U^i) = ( 27r / a )5 \Ui\ 2 d t 1 , where 
6\ is the phase of Ui, and the higher order terms with \m\ > 1 have been dropped. Therefore, 
J(t) is proportional to the absolute square of the potential and to the rate of phase change 
in the potential. 

To summarize, we have considered a ratchet model in which an overdamped particle 
is adiabatically driven by an asymmetric potential which varies periodically both in space 
and in time. The problem is treated by transforming the time-dependent Fokker-Planck 
operator to a time-dependent hermitian operator. The spatial periodicity validates the use 
of adiabatic theorem, which allows us to develop an adiabatic perturbation theory based on 
the instantaneous eigenfunctions of the hermitian operator. There are two limiting regimes, 
determined only by which of the potential energy and the thermal energy dominates. For 
each of these limiting regimes, we have discussed the adiabatic condition, for which the 
diffusion constant D sets the time scale. We then calculate the ensemble average of the 
particle velocity, which is independent of D. The necessity of parity symmetry breaking 
is reflected in that otherwise the relevant integrals vanish. An analytical expression of the 
average velocity is obtained in terms of the Berry curvature and Berry phase of the Bloch 
states. In the deterministic or tight-binding limit, i.e. when the potential dominates the 
thermal energy, the time average of the average velocity is found to be quantized as a Chern 
number, with exponentially small corrections. This analytical result confirms and extends 
a previous numerical observation, and discloses the analogy with the problem of quantum 
adiabatic transport. On the other hand it provides, in a classical statistical system, a 
physical example of the first Chern characteristic class, which has been important in gauge 
field theory [[13] and in quantum condensed matter physics P HT0] , |T^| . We have also discussed 
the kinetic limit, in which the thermal energy dominates the potential. Throughout our work, 
we have emphasized the effects of geometry and phases. 

We are indebted to Ben Simons for suggesting the analogy between ratchet and quantum 
adiabatic transport, as well as useful discussions. We also thank David Thouless and Dima 
Khemelnitskii for useful discussions. Y.S. thanks Peter Littlewood for hospitality. Q.N. was 
supported by the Welch Foundation, the Texas Advanced Research Program, and the US 
Department of Energy under Grant No. DE-FG03-02ER45958. 
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FIGURES 




FIG. 1. An illustration of a ratchet potential V(x,t) which is periodic, but asymmetric, in 
space. In our consideration, it is also periodic in time. 
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